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Abstract 



We propose a procedure based on phase equivalent chains of Darboux trans- 
formations to generate local potentials satisfying the radial Schrodinger equa- 
CD ■ tion and sharing the same scattering data. For potentials related by a chain 

\ of transformations an analytic expression is derived for the Jost function. 

I It is shown how the same system of S'-matrix poles can be differently dis- 

p • tributed between poles and zeros of a Jost function which corresponds to 

different potentials with equal phase shifts. The concept of shallow and deep 
^ \ phase equivalent potentials is analyzed in connection with distinct distribu- 

^ \ tions of poles. It is shown that phase equivalent chains do not violate the 

C7^' Levinson theorem. The method is applied to derive a shallow and a family 

of deep phase equivalent potentials describing the ^5*0 partial wave of the 
^ , nucleon-nucleon scattering. 

Typeset using REVT^ 

03.65.Nk, 11.30.Pb, 21.45.+V, 13.75. Cs 

I. INTRODUCTION 

The construction of potentials sharing the same scattering data, called in the literature 
phase equivalent potentials, is one of the most interesting applications of supersymmetric 
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(SUSY) quantum mechanics [|T| to the inverse scattering problem (see e.g. As it 

is known, the SUSY approach, when restricted to the derivation of new exactly solvable 
quantum problems, is basically equivalent to the Darboux transformations method (see e.g. 
0]). Therefore we use SUSY and Darboux transformations as synonyms. 

To our knowledge, the simplest way of constructing phase equivalent potentials by using 
Darboux transformations, is due to Sukumar [^. It consists in a successive application of 
two transformations with equal factorization constants. The method has been considerably 
developed, up to an arbitrary modification of the discrete spectrum [^]. It has been extended 
to coupled channels |0] and used to derive deep and shallow phase equivalent potentials de- 
scribing nucleon-nucleon and nucleus-nucleus PJTO| collisions. In these applications both 
deep and shallow potentials have been constructed such as to share the same scattering data 
at all energies. Shallow potentials possess only "physical" states while the deep potentials, 
in addition to the physical states also have "unphysical" bound states which simulate the 
effect of the Pauli principle (for details see e.g. |TT|-P!^). 

Based on Darboux transformations, here we propose a new way to derive deep and 
shallow phase equivalent potentials. Our aim is to improve over the practical procedures 
of Refs. and |[^-[T0[1. By replacing a chain of Darboux transformations by an equivalent 
A^th order transformation and using a Wronskian formula recently obtained in (see also 
||T5| ) we use a determinant approach to construct phase equivalent potentials. When the 
transformation functions are sufficiently simple (e.g. they are solutions of the free particle 
Schrodinger equation) this approach happens to be more efficient than the one of previous 
works |15|-[10|. The advantages of our method will be clearly seen in its application to the 



nucleon-nucleon scattering potentials considered here. In particular the potential tail will 
be cured from the unwanted oscillations obtained in Ref. 0. 

As it is known (see e.g. ^) there are both deep and shallow potentials corresponding to 
the same value of the partial wave /. The value of / defines the long distance behaviour of 
the potential 0. The short distance behaviour 

V(a:)~^±±ll, .^0 (1.1) 

depends on the parameter u which we shall call the strength of the singularity at short 
distances or simply the singularity strength. Following Ref. we generally assume that u is 
different from /. In this work we consider only l-preserving transformations and restrict our 
discussion to the / = case, although most of the derived formulae are valid in the general 
case. In a subsequent paper |T6[ we shall consider l-changing transformations which will 



allow to introduce higher partial waves. 

We assume that no Coulomb interaction is present in the potential i. e. we are dealing 
with the Schrodinger equation on the half G [0,oo), with the potential V{x) sat- 

isfying the condition ( |1.1| ) and at long distances going to zero faster than For every 
eigenvalue E the Schrodinger equation with such potentials has two linearly independent 
solutions with different behavior at short distances, one which behaves as ~ x'^, called reg- 
ular solution^ and the other as ~ x"*^"^, called singular (or irregular) solution. In general 
Darboux transformations, which use irregular solutions as transformation functions, do not 
preserve the property of a solution to remain regular at the origin. As a consequence, the 
change in the Jost function and the phase shift becomes difficult to control since one has 
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to know both the Jost solution of the initial equation and the derivative of this solution at 
the origin. As far as we know this problem has not been yet solved in the general context. 
But we have found that there exist particular chains of transformations where the derivative 
of the Jost solution of the initial equation does not appear in the expression of the Jost 
function for the transformed Hamiltonian. In this case one can get a closed expression for 
the transformed Jost function in terms of the initial Jost function and a rational function 
of the momentum k only. Potentials having Jost functions related in such a way are known 
in the literature as potentials of the Bargmann class (see e.g. 

Once we know how the Jost function is changed by a chain of Darboux transformations 
we can select chains that produce potentials with the same ^-matrix and hence with the same 
phase shift. Such potentials are called in the literature phase equivalent although a more 
appropriate term could be isophase potentials since they give exactly the same phase shift at 
all energies. Here we shall use both terms. So, to get phase equivalent potentials we propose 
to use phase equivalent chains instead of an integral form of Darboux transformations used 
in |5|-[T0|. We start from the same exactly solvable reference potential and the action of each 
chain will result in potentials always having the same S'-matrix. Finally we note that this 
approach is basically equivalent to the inverse scattering method where the kernel of the 
Fredholm operator involved in the Gelfand-Levitan-Marchenko equation is degenerate or 
to a supersymmetric transformation [^-10|, but it is more direct and has more computational 



advantages. 

The paper is organized as follows. In the next section we recall the Darboux transforma- 
tion method in the context of scattering theory. In Sec. Ill we introduce special chains for 
which a closed expression for the transformed Jost function is derived. Then we select chains 
which lead to families of phase equivalent (or isophase) potentials. Sec. IV is devoted to an 
application of our method to the derivation of deep and shallow phase equivalent potentials 
for the nucleon-nucleon scattering. Conclusions are drawn in the last section. Details of 
important analytic calculations are gathered in the Appendix. 

II. DARBOUX TRANSFORMATION METHOD 



The method of Darboux transformations, well known in soliton theory [|T7|, finds now 



more and more applications in other fields of theoretical and mathematical physics [jTS 
In this section we recall the definition of Darboux transformations and give their main 
properties necessary in subsequent sections. 

Essentially, the method consists in getting solutions (f of one Schrodinger equation 
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d 

hi^ = E^, hi = -— + Vi{x), (2.1) 



when solutions of another equation 



ho^ = Eij, ho = -— + Vo{x) , (2.2) 



are known. This is achieved by acting with a differential operator L of the form 

ip = Lip, L = —d/dx + w{x) , (2.3) 
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where the real function w{x), called superpotential, is defined as the logarithmic derivative 
of a known solution of (|2.2| ) denoted by u in the following. One has 



w = u (x) / u{x) , hou = au , (2.4) 

with a < Eq, where Eq is the ground state energy of if it has a discrete spectrum or the 
lower bound of the continuum spectrum otherwise. The function u is called transformation or 
factorization function and a its factorization constant or factorization energy. The potential 
Vi is defined in terms of the superpotential w as 

Vi{x) =Vo{x) -2w'{x) . (2.5) 

The operators L and = d/dx + w{x) establish a one-to-one correspondence between 
two dimensional subspaces of solutions of Eqs. ( |2.1| ) and ( |2.2| ) with a given E ^ a. When 



E = a the transformed equation (|2.1| ) has v = 1/u and v as linearly independent solutions. 



The latter can be found by using the Wronskian property W{v, v) = 1 which gives 

v~^{x)dx , (2.6) 



Xo 



where xo is an arbitrary constant. Hence, the knowledge of all solutions of the initial equation 
provides the knowledge of all solutions of the transformed equation and in particular of all 
physical solutions in the usual quantum mechanical sense. 

Since the above procedure does not depend on a particular choice of the potential Vq the 
transformed Hamiltonian hi can play the role of the initial Hamiltonian for the next transfor- 
mation step. In such a way one gets a chain of exactly solvable Hamiltonians ho, hi, ... hjsi 
with the potentials Vq, Vi, ... Vn- To avoid any confusion we mention that everywhere, 
except for especially mentioned cases, we shall use subscripts to distinguish between quanti- 
ties related to different Hamiltonians, h^, hi,. . . and shall omit them when discussing general 
properties regarding all Hamiltonians. 

A first order Darboux transformation operator Ljj^i as defined by (|2.3|) relates two 
Hamiltonians hj and hj+i. If one is not interested in the intermediate Hamiltonians 
hi, . . . , h^-i and all factorization energies are chosen to be different from each other the 
whole chain may be replaced by a single transformation given by an A^th order transforma- 
tion operator denoted by L^^\ defined as a succession of first order transformations. A 
compact representation of this operator is given by ||T^ 

ilj^i^x, k) = L^^'^M^, k) = W{ui, ...,un, M^, k)) W-\ui, . . . , Mat) , (2.7) 

where iIjo{x, k) is a solution of Eq. (|2.2| ) corresponding to the energy E = k"^ and ipN^x, k) 
satisfies 

hMtpN{x,k) = EtpN{x,k) , E = k'^. (2.8) 

The transformation functions, although labeled by a subscript, are eigenfunctions of the 
initial Hamiltonian 
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hoUj{x) = ajUj{x) , ttj = — a| . (2.9) 

They should be chosen such as the Wronskian W{ui, . . . ,un) is nodeless in the interval 
(0, oo). This condition guarantees the absence of singularities in the potential 

VN = Vo-2-f-logW{u,,...,UN) , (2.10) 



defining the Hamiltonian /ijv of (|2.8| ), inside this interval. The formula ( |2.7| ) is valid for any 
E = k'^ except for k = iaj (j = 1, . . . N). For these values of k the corresponding solutions 



are 



^^{x,zaj) = W^^\u,,...,UN)W-\u,,...,UM), j = l,...,N, (2.11) 

where W^^\ui, . . . , mat) is the (A^ — l)th order Wronskian constructed from ui, . . . ,un except 
for Uj. 

In general, a chain of transformations can have a number of elements corresponding 
to transformation functions with identical factorization constants. In the simplest N = 2 
case where the first transformation is implemented by the function ui = u and the second 
by f = 1/m one gets back the initial potential Vq. This means that the operator L"*" realizes 
the transformation in the opposite direction. More generally, the second transformation can 
be made with the linear combination vi = cv + v where v is given by (|2.6|) . The action of 
such a chain leads to an eigenfunction of /12 of the form 

X 

u( x] f 

ij2{x, k) = ^o{x, k) - — / u{t)Mt^ k)dt, (2.12) 



where 



xo 



X 

W2{x) =c + J u^{t)dt . (2.13) 



xo 



Choosing xq = (or xq = 00) together with the condition c > (or c < 0) guarantees a 
nodeless W2{x) inside the interval (0, 00). The integral representation (|2.12| ) based on two 
subsequent transformations with the same factorization constants was used in Refs. [^-[T^l 
for getting phase equivalent potentials. We now describe an alternative approach which has 
the advantage that the above N = 2 chain of transformations can be included in a larger 
chain sub chain. 

The potential V2(x) is a particular case of (|2.10| ) where W{ui, . . . , un) is given by W2{x) 



of ( p. 131 ). It has been shown |T^ that the function W2 is nothing else but the limiting value 



of the ratio W (ui, U2) / {ai — ^2) when U2 — > Ui. This fact, together with a property of 
the Wronskians established in |T^|T^, which allows to reduce the order of the determinants 



involved in calculations, implies that one can use the same formulae (|2.7|) and ( 2.10|) also in 



the case of pairwise identical factorization constants but the Wronskians should be calculated 
with the help of another representation. If the Wronskian in ( |2.7| ) is of even order, N = 2m, 
it becomes proportional to an m x m determinant 
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(2.14) 



(2.15) 



(2.16) 



The product ]^(ai — aj) introduced above contains only the factors (a; — aj) appearing in 

Wof,, or For identical factorization constants the function Wa,^, oi has to be calculated from 
the formula ( ^^.13] ). 



III. PHASE EQUIVALENT CHAINS OF TRANSFORMATIONS 

Our method is based on a choice of transformation functions which allows to transform 
the Jost function through a special chain of Darboux transformations. We follow the usual 
definitions of the scattering theory which can be found for example in Refs. 



A. Jost function for a special chain of transformations 

Let us choose the following set of eigenfunctions of the Hamiltonian ho 

Vi{x) , . . . ,V^{x) ,U^+i{x) ,V^+i{x) ,U„,{x) ,Vn{x) , > , 

hoUj{x) = —a^Uji^x) , hoVj{x) = — 6|fj(x) , Ima^ = Im6j = 



(3.1) 
(3.2) 



as transformation functions for the Darboux transformation of order N = 2n + u. The role 
of the first u functions Vj{x) which are absent in (|3.1| ) for i/ = will be clear below. Let 
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the a's and the 6's be different from each other and choose them such as the corresponding 
factorization energies be smaller than the ground state energy of Hq when it has a discrete 
spectrum or less than the lower bound of the continuous spectrum otherwise. 

We distinguish between the functions u and v from their behaviour at the origin. The 
functions v are regular, Vj{0) = 0, and hence are uniquely defined up to a constant factor 
which is not essential for our purpose. They form a regular family and are exponentially 
increasing at large x and therefore it is natural to restrict all 6's to be positive. The functions 
Uj are irregular at the origin, Uj{0) ^ 0, and from a singular family. They are defined as 
linear combinations of two Jost solutions f{x, ±k): 

Uj{x) = Ajf{x, —iaj) + Bjf{x,iaj) . (3.3) 

When Aj ^ the functions Uj increase exponentially at large x provided aj > 0. At fixed 
Qj they form a one-parameter (Aj/Bj) family. The case Aj = 0, Bj ^ 0, with > is 
equivalent to the case Aj ^ 0, Bj = 0, with aj < 0. They correspond to a decreasing 
exponential function ( |3.3| ), uniquely defined (up to a constant factor). In this case in order 
to stress the asymptotical decrease of this function we allow aj to be negative and set, 
Bj = 0, Aj 7^ 0. So, for Uj > the functions ( |3.3| ) have an increasing asymptotic behavior 
and for aj < their asymptotic behavior is opposite. It is useful to note that if Aj/Bj = 
F^iiaj) / Fo^—iaj) the function ( p.3| ) becomes regular. In the chain of functions ( |3.1| ) we shall 
allow to move some functions from the singular family to the regular family. This preserves 
the number of functions but changes the value of v. The effect of such a modification of 
a chain will be discussed below. The particular case where the initial potential is zero and 
Aj = Bj has been considered earlier in Ref. |^ . 



In the Appendix we show that the chain of transformations with the transformation 
functions ( p.l|) increase the singularity parameter of the initial potential by u units, z/q — »• 
= ^0 + and the transformed Jost function F/v(fc) is related to the initial function Fo(/c) 

by 

^"(*) - ^o(^-) n n f^f ■ 

j = 1 j = 

Since a Jost function is analytic in the upper half of the complex fc-plane (see P,^) all 6's 
must be positive whereas the a's can have any sign, which agrees with our sign convention 
for ttj and bj and every positive aj corresponds to a discrete level Ej = —a^ of V^- A 
remarkable property of the relation (|3.4| ) is that it is independent of the ratio Aj/Bj with 
Aj and Bj from the definition ( |3.3| ). But the potential does depend on this ratio. Hence, 
keeping -F/v(/c) fixed but varying Aj/Bj such as Uj remains irregular we can obtain a family 
of isospectral potentials. Nevertheless the potential Vn belongs to this isospectral family 
only if Aj ^ and the function Uj{x) increases at large x. The choice Aj = is equivalent 
to changing the sign of aj and leads to losing the level E = —a^ in the potential V^. This 
property does not have an analogue in the usual scheme of SUSY transformations [jl8 . 



B. Phase equivalent chains 

Now we turn to the construction of phase equivalent or isophase potentials. Such poten- 
tials have identical S'-matrices which are defined in terms of the Jost function by 
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(3.5) 



Therefore if for a real k two Jost functions differ from each other by a real factor they corre- 
spond to isophase potentials. In particular, these potentials can have equal Jost functions. 
Taking into consideration the discussion at the end of the previous subsection we see that 
by changing the ratio Aj/Bj, which corresponds to a discrete level of the potential Vjv, we 
can get not only isospectral but also isophase potentials. As it is clear from the Appendix, 
the potentials are isospectral provided the function uj remains singular. If it moves from 
the singular to the regular family, the transformed potential does not belong to the family 
of isospectral potentials any more but as we show below it remains in the family of isophase 
potentials. 

Let us suppose for example that the function Ui{x) becomes regular from being singular. 
Then for a regular Ui the Jost function of the potential Vn{x) corresponds to (|3l^) with 
z/ = and for a singular ui it should be calculated by the same formula where u = 2 and 
bi is replaced by Oi. From here it follows that the ratio of these two Jost functions is equal 
to fc^/ (/c^ + a1) which is real for a real k and hence these two chains of transformations are 
isophase. It is also clear that if u functions from the singular family move to the regular 
family (i.e. u zeros of the Jost function are transformed into its poles) the potential Vn loses 
ly discrete levels, and acquires an additional singularity strength equal to 2i/, but remains 
in the family of isophase potentials. Furthermore, this procedure is completely reversible. 
If one has a singular potential with u > 2 then one of the poles of its Jost function can be 
transformed into a zero in the upper half of the complex fc-plane which leads to the creation 
of a discrete level and a decrease of the singularity strength u u — 2. This procedure can 
be repeated either up to the point where the potential loses all discrete levels and all free 
parameters, becoming uniquely defined by its phase shift (see e.g. 0,|§]) or up to the point 
where u becomes equal to zero or one. To avoid additional singularities inside the interval 
(0, oo) one has to take these poles such as the additional level be located below the existing 
discrete levels. The latter restriction is well known in the SUSY approach 



The above method of deriving isophase potentials is not unique. As it is known one 
can create a new energy level independent of the existence of poles of the Jost function, 
at any desirable position even in the middle of the continuous spectrum, and preserving 
the S'-matrix. For this purpose one can use a two-step Darboux transformation with the 
same factorization energy After a new level has been created the singularity of the 

potential loses two units u u — 2 [0. This procedure can be repeated up to a limit where 
the singularity parameter becomes equal to either zero or one. For a regular potential or a 
potential with u = 1 no possibility exists to create a level by a phase preserving procedure. 
In this context the smallest possible value is z/ = 2 when the potential with an additional 
discrete level becomes regular after transformation. According to Ref. 0] a chain of two 
transformations with coinciding factorization constants is equivalent to considering the limit 
CLjji_i — > — CLj and using the Wronskian formulae ( p.l4|) and (|2.15| ). Without loosing generality 



let us consider 02 — >■ — Oi for our system of transformation functions ( |3.1| ) with i/ = 0. In 
this limit the Jost function (13.41) where we set z/ = 



Fr,{k) = Fo{k)—^^±p—— ff . (3.6) 

^ \k + ibi){k + ib2) ^-J-^k + ibj ^ ^ 
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has no singularity. The corresponding potential Vn can be considered as a result of the appli- 
cation of the above two-step phase preserving SUSY transformation on a singular potential 
V^r_2 generated by the system of transformation functions Vi,V2,U3,V3, . . . ,Un,Vn whose Jost 
function -F/v_2 is given by the same formula (|3.4|) with u = 2. Finally we conclude that V/v-2 
is isophase with Vn since the ratio Fn/Fn-2 = {k'^ + ci\)/k'^ is real. 



C. The S"- matrix poles 

Once the Jost function F{k) is known one can calculate the phase shift 5{k) by using the 



definition ( |3.5| ). The asymptotic behavior of the scattering wave function of the potential 
Vn with I = and singularity strength u at long distances is 



ipiyix, k) ^ sm{kx — ^un + Sj^fik)), x 



oo 



(3.7) 



In our approach F]sr{k) differs from -Fo(^) by a rational function of momentum k and therefore 
the expression for S^ik) becomes rather complicated when the number of transformation 
functions is sufficiently large. An equivalent expression which is more convenient for practical 
calculations is 



^Nik) = So{k) — arctan(A;/aj) — a.Ycta.n{k/bj] 



(3.8) 



in agreement with the result obtained in From ( p. 81 ) it follows that the S*- matrix acquires 
additional poles at aj and bj through the application of the transformation L^^^ consistent 
with the discussion of the previous section. 

For / = the effective range expansion is (see e. g. [pl[] ) 



lim k cot S(k) 



1 1 ,2 

- + T^rok^ 
a 2 



Pr-^k'' + ... 



where a is the scattering length and tq effective range. Expanding the phase shift 
power series one obtains 



(3.9) 



m a 



1 



1 



lj = u+l 



(3.10) 



and 



,(JV) 



-3 



(3.11) 



These formulae can be used to calculate theoretical values to be compared with experimental 
results, as for example those obtained in the neutron-proton scattering experiments. But 
first one has to fit the parameters aj and bj in (|3.8|) such as to reproduce the phase shift 
data as close as possible. As a result of the fit one gets a sequence of n+ positive and 
negative poles of the S'-matrix all having absolute values below the pole associated with the 
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ground state of ho if this has a discrete spectrum or below the lower bound of its continuum 
spectrum otherwise. (We shall comment below on the cases when this is not so.) Every 
negative pole is associated with a singular transformation function decreasing at infinity 
and belonging to the singular family of (|3.1| ). This is just the Jost solution of the initial 
equation at kj = iaj which is uniquely defined. A positive pole corresponds to an increasing 
transformation function (|3.3|) , which can be either regular or singular at x = 0. When 
regular it is uniquely defined, whereas singular functions form a one-parameter family. 

According to the choice of transformation functions (|3.1| ) our method is restricted by the 
condition ?7,_ < n+ and one has to distribute the functions related to positive poles among 
regular and singular families such as the number of singular transformation functions does 
not exceed the number of regular ones. This means that at least n_ functions related 
with positive poles should be regular. For the remaining positive poles there are different 
possibilities either related to the regular {vj} family or to the singular {uj} family. All these 
possibilities correspond to the same phase shift and hence will give different phase equivalent 
chains of transformations. When all the remaining positive poles are associated with the 
regular family the potential will acquire the largest singularity strength z/ = — r2_ and 
according to the formula (p.4|) will have no additional discrete levels with respect to Vq. In 
this case every function from the system ( |3.1| ) is uniquely defined and the potential has 
no free parameters. This agrees with the well-known statement of the inverse scattering 
method (see e.g. @J3) that in the absence of discrete spectrum a local potential is uniquely 
defined by its phase shift. Such potentials are called in the literature shallow potentials. 

When a maximal number of positive poles (equal to 7 = (n+ — nJ)/2 for even and 
7 = (n+ — n_ — 1)/2 for A^ odd) is included in the singular family every such pole corresponds 
to an additional discrete level of Va? and this potential will acquire the smallest singularity 
strength v = Q (regular potential) for A^ even and v = 1 for A^ odd and 7 additional discrete 
levels. Every positive pole aj associated with a singular function Uj will correspond to a 
discrete level Ej = —a^ of Vn- According to ( ^.31 ) the function uj has a free parameter Aj/ Bj . 
Thus we obtain a 7-parameter family of isophase and isospectral potentials. Such potentials 
are called in the literature deep potentials. Of course there are also intermediate possibilities. 
When = ?t,_ or n_|_ = n_ + 1 one can construct only a deep (regular if n_|_ = n_) potential 
without additional discrete levels. For this potential phase preserving transformations which 
at the same time do not change the angular momentum / are impossible. 

When 1/ = n_|_ — n_ > 2 one can construct a 27-parameter family of isophase potentials 
by introducing 7 energy levels into the shallow potential. For this purpose one has to use a 
chain of A^ + 27 transformations, phase equivalent to the initial chain of A^ transformations, 
to which 27 transformation functions with pairwise equal factorization constants = 6| 
should be added. To calculate the Wronskian of such a system one has to use either the 
formula ( p.l4| ) or (p.l5|) . The potentials form a 27-parameter family where 7 energy 
levels play the role of parameters. The other 7 parameters correspond to the constant c of 

It may happen that the absolute values of some poles is above the position defined 
by an excited state energy of Hq. Then using a phase preserving procedure, we first have 
to delete all such levels including ground state level from the spectrum of Hq and then 
use the above described technique. This means to include additional pairs of functions 
with equal factorization constants in the system of transformation functions (p.ll). The 
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corresponding energies have to coincide with the positions of the discrete levels to be deleted, 
which evidently does not change the fitted phase shift. The Wronskian of the entire system 
of transformation functions has to be calculated by using either ( p.l4| ) or ( p.l5| ). 



When the number of negative poles exceeds the number of positive ones, the problem 
becomes more sophisticated. In some cases it can be solved by introducing transformations 
changing the orbital angular momentum I. We are planning to discuss these problems in a 
subsequent paper []TB|. 



In this subsection all potentials have an exponential decrease at infinity, provided the 
initial potential corresponds to / = 0. Therefore all transformed potentials also correspond 
to / = 0. An useful remark is that when fitting an experimental phase shift one can consider 
either real poles of the S'-matrix or pairwise mutually conjugated poles. In either case 
they should be associated with the regular family of the system The corresponding 

Wronskian is either real or purely imaginary, since the transformation functions are mutually 
conjugated. The potential difference Vjy — Vq resulting from ( |2.1(J| ) is real in either case. 

D. The Levinson theorem 



Here we show that the formula (|3.8| ) is in agreement with the Levinson theorem in 
its generalized form |^ (see also P). For a potential satisfying the condition ( |1.1| ) and 
decreasing faster than at long distances this theorem states that 

(5(0)-5(oo) = [A/'+i(z/-/)]7r , (3.12) 

where A/" is the number of bound states I the partial wave and u the singularity strength. In 
our case we have to set / = 0. Suppose that this theorem holds for the initial potential Vq. 
We shall show that in this case it also holds for the transformed potential Vat. Since the right 
hand side of ( |3.12| ) is a linear form in the variables Af and u it remains to show that this 
formula is valid for the case where 5(0), ^{oo), M and v represent differences between the 
same quantities associated with the potentials Vq and Vat. In what follows in this subsection, 
by these quantities we understand the corresponding increments. 

The potential Vat is supposed to be obtained from Vq by an /-preserving chain of trans- 
formations with transformation functions ( p.l| ). When all additional positive poles of the 
S'-matrix are related to the regular family, the potential Vat has the same number of bound 
states as Vq but its singularity strength is z/at = = r;,+ — n_ . Using ( |3.8| ) one finds that 
5(0) = 0, (5(cxd) = — 7r(n+ — nJ)/2 which is in the agreement with ( p.l2|) at A/" = 0. If 
now we transform M < (n+ — nJ)/2 regular functions into singular ones the potential Vn 
acquires M additional discrete levels and loses 2N units of singularity strength which also 
is in agreement with the relation ( 3.12| ). Hence, we conclude that /-preserving chains of 



transformations with transformation functions ( p.l|) do not violate the Levinson theorem. 



IV. APPLICATION TO NEUTRON-PROTON ELASTIC SCATTERING 

Here we apply our method to derive shallow and deep phase equivalent potentials, both 
reproducing the neutron-proton {np) experimental elastic scattering phase shift for the 
^5*0 partial wave in the laboratory energy interval 0-350 MeV. The experimental values 
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taken from Ref. p3| are displayed in Table I. Our choice Vo{x) = which allows to use 



solutions of the free particle Schrodinger equation simplifies the calculations considerably. 
The theoretical phase shift also displayed in Table I under 6q is obtained from the formula 
( p.8|) with 6(){k) = and u = where we had to fix n optimalUy and search for aj and bj 



(j = 1,2, ...n). By using the standard Mathematica package NonlinearRegress we found 
that we need six S-matrix poles in order to fit the experimental phase shift very closely. 
Thus we have N = 2n = 6 and the poles are 



ai = -0.0401, as = -0.7540, as = 4.1650, 
bi = 0.6152, 62 = 2.0424, 63 = 4.6000, 



(4.1) 



in fm ^ units. In the notation of Subsection III C we have n_ = 2 and n+ = 4. From 
Eqs. (|3l0|) and (|3Jl|) with l/a^") = and r^^^ = we obtain a^^'^ = -23.7032 fm and 



r*-^-* = 2.6235 fm. The comparison with the recommended values aexp = —23.721 ± 0.020 



fm and r^xp = 2.658 ± 0.062 fm shows that the phase shift is sufficiently well described 
by the system of poles ( |4.1| ). The negative poles correspond to transformation functions of 
type Uj = exp(ajx). Positive poles can correspond either to regular solutions Uj = sinh(aja:) 
or to 



u 



Aj exp{ajx) + Bj exp{—ajx) , (4.2) 



as implied by (|3.3|) . According to the previous section the poles bj correspond to regular 
solutions Vj = sinh(6jx). 

We can now derive isophase potentials associated with Sq by using ( p.lQ ). The system 
of poles ( |0| ) can correspond either to a one level potential (deep potential) with Eq = —al 
or to a shallow potential which has no discrete levels. The latter is the desired case and can 
be obtained when the function M3 is a sinh function. Fig. 1 shows our shallow potential Ve 



obtained in this way together with the Reid68 potential 26 



Vneidesix) = [-10.463 exp(-0.7x) - 1650.6 exp(-2.8a;) + 6484.2 exp(-4.9x)]/(0.7a;) , (4.3) 



and the potential of Ref. The dots represent an updated version |^3[ of the Reid soft core 
potential ||26|. In contrast with Ref. our potential has a correct short range behaviour 
Vq{x) 6x~^ corresponding to the case u = rij^—n- = 2 (the potential of Ref. has u = 1). 
In Fig. 2 we present the tail of the absolute value of various potentials in a logarithmic scale. 
One can see that Vq has a smooth behaviour consistent with the one-meson exchange theory 



as it is the Reid's potential ||2^, in contradistinction to the SUSY potential from Fig. 3 of 



Ref. IP which displays an undesired oscillation around r = 2 fmQ. Interestingly, our potential 
is very close to the Reid68 potential ([4.3|) in a large interval up to 5 fm. 

Actually the choice (|4.2| ) for gives an isospectral and isophase family of Vq potentials 
with the bound state energy Eq fixed. In Fig. 3 we display a few deep potentials, members 
of this family and isophase with the shallow Vq potential. These potentials differ from each 



^The potential drawn in Figs. 1 and 2 under Ref. |^ represents the output of Eq. (32) of this 
reference. 
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other only by the value of the coefficient ^3 = ^43/ B3 in ( [4 .21) . One can see that by increasing 
A'^ the potential well shifts to the right without changing its depth. The bound state with 
£"0 = — ag follows the movement, being localized at larger and larger values of x. A similar 
behavior has been observed by Bargmann already in 1949 (see also p8|). When A'^ 



approaches —1 the well concentrates near the origin and becomes very narrow and deep. 
The ground state becomes well localized near the origin. At very short distances a barrier 
appears before the well which grows when — > — 1. In the limit A'^ = —1 the well becomes 
infinitely narrow and icompletely collapses, being replaced by an infinitely high barrier which 
becomes a repulsive core, as for a shallow potential. 

One can also construct a deep potential isophase with the shallow potential Vq but having 
a discrete level at a desired place. First one has to delete the bound state Eq = — Og of Vq 
following the procedure described in Sec III. Next one has to use a 2-parameter family of 
isophase potentials (Subsection III C). This is possible by adding two new transformation 
functions, which gives a total of 8 functions. The new functions M4 and f 4 must have 04 = —64 
(a4 < 0) in order to give Ss{k) = Sei^k). It is enough to take M4 = exp(64x) and calculate 
the Wronskian of the added pair from ( |2.13| ). When |a4| = 03 the potential behaves like 
a deep Vq potential with the energy E = — Og. By varying 04 one can change the depth 
of the well. The family of Vg potentials has in fact two parameters: the energy E = —a1 
and the parameter c of ( |2.13|) . One can select a potential from this family by fixing these 
two parameters. Here we look for a deep potential which describes the nucleon-nucleon 
interaction according to microscopic theories based on the quark structure of nucleons (see 
e. g. Refs. [0 or p9|)- wish to obtain a potential close to the deep potential Vk of Ref. 
29| , given by the formula 



Vk = -1106.21 exp(-1.6x2) - 10.464 exp(-0.7x)/(0.7x) [1 - exp(-3x)] , (4.4) 
where x is given in fm. We varied 04 and c such as to be as close as possible to ([4.4[). 



From a least square fit we obtained c = —0.155 fm and = —3.7944 fm~^ which gives E 
= 596.42 MeV. In this way we obtained the potential Vg which is compared to Vk of ( |4.4|) 
in Fig. 4. One can see that the two potentials are quite close to each other. The bound 
state of Vk is 442. MeV which is not much different from our value. The form of Vk was 
chosen on qualitative grounds related to the approximations used in microscopic studies 
of the nucleon-nucleon interaction [|12|. The potential Vg of Fig. 4 is as legitimate as Vk 



inasmuch as it correctly describes the experimental phase shift. The difference between the 
two potentials is partly due the difference in the fitted phase shifts and the Eiab interval 
chosen for the fit in each case. If instead of fitting the phase shift to the experimental values 
we were fitting it to the phase shift corresponding to the potential (f4.4|) , the difference would 
have been smaller. 



V. CONCLUSIONS 



In this study we propose an alternative method to derive families of phase equivalent 
(isophase) local potentials based on the application of an iVth order Darboux transformation 
or an equivalent chain of ffist order transformations. In practical applications this method 
is simpler and more efficient than the conventional SUSY approach of Refs. p|-pi]|]. First 
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one has to look for imaginary poles of the S-matrix, i.e. one has to find the parameters Oj 
and bj of the phase shift S^lk) of ( p.8|) which fit a particular experimental phase shift. One 
can start with the lowest order = 2 in ( |3.8| ) and increase it until a good fit of the phase 
shift is achieved. The second step consists in finding an appropriate mapping between the 
S-matrix poles found in this way and the system of functions which give a nodeless 
Wronskian. In terms of these functions the formula (|2.10| ) gives phase equivalent potentials, 
among which one is shallow and the others form families of deep potentials from which one 
can select a particular one for a specific physical problem. We have applied our method to 
derive potentials describing the neutron-proton experimental phase shift for the ^5*0 partial 
wave. We obtained a shallow potential very close to the Reid's soft core potential and its 
deep phase equivalent partner is close to potentials resulting from microscopic derivations 
of the nucleon-nucleon interaction. In order to obtain phase equivalent potentials describing 
higher partial waves /-changing Darboux chains of transformations are necessary. This is 
the scope of a subsequent study. 



APPENDIX A: 

Here we prove the formula ( [3.41 ). The proof is split into two steps: 1) Taking an initial 
potential with i> = uq = we prove ( p.4| ) for this case. 2) We introduce a number of regular 
solutions to the previous case which corresponds to z/ > 0. Finally we comment on singular 
initial potentials. 

Let us consider a chain of N = 2n Darboux transformations for a potential Vo{x) with 



v = = and transformation functions given by (^) with z/ = 0. Since all a's and fe's 
are different from each other this chain is equivalent to an A^th order transformation. The 
transformed potential Vn is given by ( p.lO| ) and the eigenf unctions ipN by ( p77| ) or ( p. 111 ). 



One can also find the Jost solution of ftTby applying ^he same operator ii^to the jSt 

solution of ho. The result is the function /Ar(x, k) = L^^^ fo{x, k), which like i/jn from 



is proportional to the ratio of two determinants. The long-distance asymptotic value of 
fiy{x,k) is proportional to the asymptotic value of the Jost solution. Therefore using the 
usual definition of the Jost function as the value of the Jost solution at the origin one 
has 

= — km — (Ai) 

[e *''^/Ar(x, /c)]^_oo 



The definition ( p. 71) and our special choice of transformation functions ( p.lD implies that 

fN{x, k) = W{Ui, Vi,..., Un, Vn, fo{x, k)) W'^{ui, Vi,...,Un, Vn) = Wf W^^ (A2) 

Therefore both for x = and x ^ oo we have to calculate the ratio of two determinants of 
orders + 1 and respectively, but of a similar structure. 

First let us consider the function /Ar(x, k) ai x = 0. Using the definition /o(0, k) = Fo{k) 
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the determinant appearing in the numerator of ( |A2| ) becomes 



Wf 



Ci 





C2 





■ ■ ■ Cn 





Fo{k) 


Cl 


Di 




D2 


■ ■ ■ 


Dn 


F{k) 


Cial 





C2CI2 





■ ■ ■ Cntt^ 





Fo{k){ikf 
F{k){ikf 


Cial 




C2a\ 


D2hl 


■ ■ ■ CfiOjj^ 


Dnhl 








C2al^-^ 










Foik)itky^~' 






C2a'r' 


D2hl^-^ . 






F{k){ikf''-^ 

















Fo{k)iikf" 



Ci...CnDi...DnFoik) 



al 



2n-2 „2n-2 



"1 
2n 



On 



,2n 



2n-2 /■i\2n-2 



,2n 



{zky 
[ikf'' 



1 

hi . 


. 1 

■ bl 




t2n-2 
■ n 



-Foik) lliC.D,) n (a? - - bl) n(«? + 

j =1 m>j j =1 



(A3) 



Here Cj = Uj{0), cj = u'j{0), Dj = Vj{0), F{k) = fo{0,k). Note that the derivative of the 
Jost function, F{k), does not appear in the final expression for Wf. The above result has 
been obtained by first using the following determinant property 



pxp 







M, 



qxq 



det Mpxp det Mq^q 



(A4) 

from the middle 



and next by expanding the two Vandermonde determinants (see e.g. 
line. 

The determinant Wn which appears in the denominator of ( [A^ ) is of even order. It is 
obtained from the above determinant by the removal of the last row and column. Its final 
form is similar to the one above. It has opposite sign and the factors depending on k are 
absent. Finally one finds 



fN{0,k) 



-irF,ik)i[ik' + 



(A5) 



Let us now consider the denominator of ([Al|), i. e. fN{x,k) for x —> 00. In this limit 
all functions entering the Wronskians of (|A2| ) tend to exponentials which cancel out with 
each other except for exp{ikx) which cancels out with the exponential present in ( [All) . 
The remaining factor is a ratio of two Vandermonde determinants of orders iV + 1 and 

respectively containing powers of aj and bj. The determinant in the numerator differs 
from that in the denominator by the presence of a column containing powers of ik and an 
additional row composed of (A^ + l)th power of ai, 61, . . . a„, 6„, ik. Using once more the 
expression of Vandermonde's determinant one obtains 
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fN{x, k)U^^ -^^^Yl{k + ia,){k + ih^) (A6) 
i=i 

Introducing ( |A5|) and (^^) in ([Al|) one gets the following relation between the initial 
and final Jost functions 

" h — ' 

f»w = F„(*)nFnr- m 

j=i J 

Moreover, since the Wronskian W^q is different from zero at x = the singularity strength 
parameter of Va? is zero and it is regular at the origin. 

Consider now the case where the first z/ functions of (|3.1|) belong to the regular family 
and calculate the singularity strength of Vn- Note that the system can be divided 

into two subsystems, one containing an equal number of singular and regular functions, and 
the other composed of regular functions only. Accordingly the A^th order transformation 
can be decomposed into a superposition of two transformations. The first transformation 
of order N — v which is based on functions paired by their regular property will not affect 
the regularity of the potential at the origin, whereas the second one of order v composed of 
regular functions only will increase the singularity parameter of the potential by one unit for 
each transformation function (see e.g. 0). So, the singularity strength parameter of Vn{x) 
is equal to z/. 

In the following we have to use the definition of the Jost function 

F(fc) = lim f^^'^ (A8) 
t^lzu — 1)!! 

We give detailed calculations for the particular case u = 2 from which a general formula 
can easily be derived. In this case both Wj = W{vi,V2,U3,V3, . . . ,Un,Vn, fo) and = 
W{vi,V2, Ms, f3, . . . , Un, Vn) are equal to zero at x = 0. Since /o(0, k) only one function 
(either vi or V2) in Wf is an unpaired regular solution and hence the zero of Wf{0) is of 
the first degree. For i> = 2 the second logarithmic derivative of Wi\i{x) at x = is of order 
z/(z/ + 1) = 6. We conclude that the zero of Wn{0) is of the third degree. To evaluate the 
indeterminate form defining the Jost function (see ( |A8| ) for u = 2 and f{x, k) = fN{x,k)) one 
can use L'Hospital's rule to get {kxf Wj/Wn = k^ Wf/x * x^Wn 3! A;^ Wj{0)/W^{0) 
when X — 0, because the first two derivatives of Wn{0) are zero. Since the derivative of an 
A^th order Wronskian W{fi, . . . , /jv) is a determinant obtained from the same Wronskian 
by the replacement of the last row composed of the derivatives fj^~^^ by the derivatives 

/j^'* one can use the property (|A4|) to factorize Vrj(O) into a product of two Vandermonde 
determinants Wyi and Wv2 of orders n and n + 1 respectively. When calculating the third 
derivative of Wn at x = one can see that it is equal to a sum of two equal determinants 
which can be factorized according to ( |A^ ) in a product of two Vandermonde determinants 
Wv3 and Wva of orders n — 1 and n + 1 respectively. We do not give here the explicit form 
of these determinants. The factor 2 coming from two equal Vandermonde determinants will 
cancel out with the corresponding factor coming from 3! and the factor 3 from this factorial 
will cancel out with (2z/ — 1)!! for = 2, the latter being present in the definition of the 
Jost function (^. Hence, {kx^Wj/WN k'^WviWv2/ {Wv2.Wv4)- When evaluating the 
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corresponding Vandermonde determinants one can see that the factors independent of k are 
exactly the same both in the numerator and the denominator. So, one gets the following 
Jost function for V^^x) 



n 



J — -J 



It is clear from here that every additional regular function in the system (|3.1|) with the 
parameter bj produces a factor k/{k + ibj) in the Jost function which corresponds exactly 
to the formula ( |3.4| ). Moreover, a regular at the origin potential admits two kinds of trans- 
formations, the ones which do not affect its singularity {u = in ( |3.1|) ) and others which 
increase its singularity {u > 0). Actually if the initial potential has the singularity param- 
eter z/q 7^ 0, three types of transformations are possible. In addition to two the previous 
types of transformations existing for the case 1^0 = transformations which decrease the 
singularity parameter are also possible. All these possibilities become clear if one takes into 
consideration the fact that all transformations are invertible (this means that from an inter- 
mediate element of a chain one can move in both directions) and an intermediate member 
of a previously considered chain is just a potential with a singularity strength different from 
zero. 
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TABLES 

TABLE L The experimental phase shift from Ref. together with the theoretical value 



6q. Eiab is related to the momentum kifm "'^) by A; = \ — " where we took m„ = 940 MeV 

V 2?i^c^ 

and he = 197.33 MeV fm. 



Eiab (MeV) 


5exp{deg) 




0.100 


38.43000 


38.422 


4.000 


64.28440 


64.362 


9.000 


60.72287 


60.717 


17.000 


55.34120 


55.253 


28.000 


49.39220 


49.291 


42.000 


43.32663 


43.270 


56.000 


38.29709 


38.291 


70.000 


33.95623 


33.988 


84.000 


30.11287 


30.167 


98.000 


26.64896 


26.711 


112.000 


23.48572 


23.546 


126.000 


20.56759 


20.620 


140.000 


17.85371 


17.893 


154.000 


15.31309 


15.337 


168.000 


12.92158 


12.930 


182.000 


10.66001 


10.654 


196.000 


8.51284 


8.4939 


210.000 


6.46736 


6.4377 


224.000 


4.51295 


4.4752 


238.000 


2.64074 


2.5979 


252.000 


.84314 


0.7984 


266.000 


-.88631 


-0.9299 


280.000 


-2.55320 


-2.5924 


294.000 


-4.16240 


-4.1942 


308.000 


-5.71816 


-5.7397 


322.000 


-7.22423 


-7.2326 


336.000 


-8.68393 


-8.6765 


350.000 


-10.10022 


-10.074 
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FIGURES 




r(fm) 

FIG. 1. Comparison between the shallow potential Vg (solid line), the Reid68 potential |26 
(dotted line), the potential of Ref. |^ (dashed line) and the Reid93 potential [^] (dots) represented 
as a function of x = r. 
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FIG. 2. The asymptotic behaviour of the absolute value of (solid line), the Reid68 potential 
(dotted line) and the potential of Ref. [0] (dashed line) in (natural) logarithmic scale. 
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FIG. 3. Family of isospectral potentials isophase with Vq. Solid line - A'^ = 0; Dashed line - 
^3 = 10^ dash-dotted line ~ A'^ = -0.95. For the latter case appearance of a barrier is clearly 
visible. 
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FIG. 4. Comparison between two deep potentials: Vs{x) (solid line) and the potential of Ref. 
(dashed line). 
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